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Abstract 

The U(l,l) and U(2) transformations realized by three-mode interac- 
tion in the respective parametric approximations are studied in conditional 
measurement, and the corresponding non-unitary transformation operators 
are derived. As an application, the preparation of single-mode quantum 
states using an optical feedback loop is discussed, with special emphasis 
of Fock state preparation. For that example, the influence of non-perfect 
detection and feedback is also considered. 

PACS numbers: 42.50.Ct, 42.50.Dv, 03.65.Bz 
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Conditional measurement offers a promising way to manipulate the state of a 
given quantum system. The basic idea is to entangle the state of the system un- 
der consideration with the state of an auxiliary system and to prepare the system 
in the desired state owing to the state reduction associated with an appropriate 
measurement on the auxiliary system. In what follows, we restrict our attention 
to travelling optical fields. The quantum states of two travelling modes can be 
entangled by mixing them at an appropriately chosen multiport. Possible basic 
transformations are the U(l,l) transformation as realized by a non-degenerate 
parametric amplifier and the U(2) transformation as realized by a frequency con- 
verter or a beam splitter. 

The aim of this paper is to generalize and unify previous work on conditional 
measurement at U(2) and U(l,l) couplers. This includes the description of the 
quantum-state transformation in terms of a non-unitary operator as well 
as possible applications such as the generation of Fock states [Q, Schrodinger 
cat-like states and photon-subtracted or photon-added Jacobi polynomial 
states or the measurement of specific overlaps |^. By combining a theo- 
retical concept for preparing single-mode quantum states by alternate coherent 
displacement and photon-adding [§] with an experimental proposal to employ 
a 2-photon down converter inside a feedback loop 0, a way is offered to pre- 
pare quantum states without nonclassical input. In particular, emphasis is being 
placed on the generation of Fock states, for which possible imperfections of the 
optical components are analysed. 

The paper is organized as follows. In Sec. |I| appropriatly factorized represen- 
tations of the unitary U(l,l) and U(2) transformation operators are introduced. 
They are used in Sec. ^ in order to derive the nonunitary transformation opera- 
tors realized by conditional measurement. In Sec. |^ the results are applied to the 
generation of specific quantum states, with special emphasis on Fock states and 
simple superpositions of Fock states. Finally, a summary and some concluding 
remarks are given in Sec. ^ 

1 U(l,l), U(2) and \J{p,q) from parametrically 
approximated three-wave-mixing 

Let us consider the transformation 

^abc = UgatcU' (1) 

of the quantum state Qabc of three travelling optical modes (denoted by a, b, and 
c), with U = e~^^* being realized by a three- wave-mixer, 

H = huJaO^a + huj^h^h + huJcC^c + /ix^^H^^^^c + c^a). (2) 

Here, t may be regarded as being the interaction time and x*^^"* corresponds to the 
second-order nonlinear susceptibility. If we assume that the mode c is prepared 
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in a coherent state, Qabc = Qab ® |7)(7|; then in the hmit x^^"* and I7I 00 
with x'^^)7 = const, the reduced density operator Q'ah = '^'^cQ'ahc reads 



Qab — UhQabUA 



A' 



(3) 



where Ua = g r"^^* is a U(l,l) transformation reahzed by a non-degenerate para- 
metric amphfier, 



Ha = huOaCi^a + huobV^h + hx^"^^ {-ia)P + 7*6a) . 
Introducing the quantities 



00 = -{^a - ^b)t, 

01 = -(7 + 7*)x^'^t, 

02 = -i(7-7*)x(%, 

where the commutation relation 



Ko = liaa^ -b^b), 
ki = \{b^a) + ah), 



(4) 



(5) 



[kj,kk] = i ^ SojkigimKr, 



l,m=0 



(6) 



is vahd = diag(l, 1, 1, —1), Sijki is the four-dimensional Levi-Civita symbol], 
and applying the respective disentanglement theorem [0, we may factorize If a 
as follows: 



where 



Ua 



g-i^^^gi</'0^0gi(¥'T+¥'fl)i^"3g2ii»i'2gi(</PT-V'fl)^3 



rp*-l ^prp\^*—a^a^—P*Ry^a^^PR*ab^prjn*Y 



(7) 



T 

p 



cosh ^^e^'''^ = cosh — hi — sinh — 
2 6 2 



sinh ?9e'^« 



+ 1^ 



sinh 



\20. 

e 2 , 



(8) 

(9) 
(10) 



and T = Te-' 2 (0 = v/0? + 0i-0^). For 0o = we have P = 1 and ^/a in (0) 
reduces to a SU(1,1) transformation operator, compare 

If, alternatively, we assume that the mode b is prepared in a coherent state, 
Qabc = ^ac ® 1/^) , then in the limit x^'^^ ~^ and |/?| — > oo with x^^"*/? = const, the 
reduced density operator g'^^ = Tr{,^^j,^ reads 



Qac = UcQacU, 



t 

c 
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where Uc = e ft"^*^* is a U(2) transformation realized by a frequency converter, 

He = TiuJaO^a + TiujcC^c + ^x^^^ (/^c^a + P*a)c) . (12) 
In order to factorize f/c, we introduce the quantities 



^1 = -{(3* + P)x^% h = |(atc + eta), 

^2 = -i(/?* - f3)x^^H, U = ^(a^c - eta) 

933 = -{ua- ujc)t, L3 = i(ata - etc). 



(13) 



where now 
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[Lj,Lk] = i^sojkiLi. (14) 



1=0 



We apply the respective disentanglement theorem [|T^ to obtain 



where 



T = cos'i^e"''^ = cos — + i— sin — , (16) 

1 Lp 1 

7c = sm'(ye^-^ = sm — , (17) 

P = e^^ (18) 

((/) = + + v^l). For (^o = we have P = 1 and ?7c in ([T5|) reduces to a 
SU(2) transformation operator, compare [|1^ . 



From (1^) and ([15|) the transformation matrices for the respective mode oper- 
ators are deduced to be 

In turn, (|TP|) and ( PUj) can themselves be used to define the four-parametric action 
of a parametric amplifier and a frequency converter. In this case we consider T, 
i?, P and T, P, V as six complex numbers that satisfy the four conditions 

|T|2 - |P|2 = |P|2 = |r|2 + |P|2 = |P|2 = 1 (21) 
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and are otherwise arbitrary. Using ( ]I9| ) and pO|), it is not difficult to verify that 

[a^a - ¥b, Ua] = [a^a + c^c, Uc] = 0. (22) 

The inverse transformations are obtained by replacing (pj with — 0j and (fj with 
-ifj in (D - (0) and (|16D - ([T|), respectively, i.e., 

UX\T,R,P) = Ua{T*,-R,P*), (23) 
Uc\T,n,V) = Uc{T\-n,V*). (24) 

and interchanging signal and idler mode leads to 

UA{b,a;T,R,P) = UA{a,b;T, R, P*) (25) 
Ucic,a;r,n,r) = Ucia,c;r*,-n*,V). (26) 

If we assume that modes b and c are simultaneously prepared in coherent 
states, Qabc = Qa ® ® |7)(7l; "then in the limit x^^^ ^ 0, \f3\ ^ oo, I7I 

00 with x^'^^Pl* = const, the reduced density operator g'^ = TibcQabc reads g'^ = 
UDQaUl, where 

= e-^^"*"^'^Z)(a) (27) 

can be written as the product of a U(l) transformation e~"^°*''^" and a coherent 
displacement D{a = iF*) = where F = Fd + F*a) and F = ix^"^^ PY^a^i^ 
— e~"^"*) (for details on the conditions under which parametric approximations 
hold, see H). 

Linear coupling of more than two modes can be reduced to a successive ap- 
phcation of two- mode couplers. For instance, as generalization of (|^) and (pOD , 
let us consider a U(p,g)-coupling of two sets of modes ai, ■ ■ -, and cip+i, ■ ■ ■, a^r 
(N = p + q), 

g-iatHap^giatHa _ ^iGH^^ ^28) 

where a is a column vector whose elements are Si, ■ ■ -, Sp, dp^^, ■ ■ ■, a]v and is a 

row vector with elements al, ■ ■ ■, a],, Sp+i, ■ ■ ■, Sat. H is a hermitian N x N matrix, 
and G is a X diagonal matrix whose upper p (lower q) diagonal elements 
are equal to 1 (—1) (see Appendix 0). Corresponding to a factorisation of e^**^^**, 
(|28D can be implemented by successive application of (^) two-mode couplers, 
each connecting two of the A^ modes by either a U(2) transformation (if the two 
modes belong to the same set) or a U(l,l) transformation (if the two modes 
belong to different sets). The special case of implementing U(A^) = U(A^, 0) by 
means of beam splitters is discussed in Note that e^*^^^** can alternatively be 



factorized into U(2) transformations and single-mode squeezing operations |T 
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2 Conditional measurement at U(l,l) and U(2) 
couplers 

Let us consider the scheme in Fig. ^ The signal mode a (index 1) prepared 



Alice 



\G) 




Figure 1: Schema of controlled quantum state engineering by conditional measurement at a 
U(l,l) or U(2) coupler realizing the transformation U = Ua oi U = C/c, respectively. If Alice' 
measurement device has detected a desired state |G), she informs Bob who opens the aperture 
S and 'stores' the pulse in the medium M until needed ||l6|] . 

in a state and the idler mode {b or c, index 2) prepared in a state \F) are 
mixed at a U(l,l) parametric amplifier or a U(2) frequency converter, and a 
device D performs some measurement on the output idler mode. (The pump 
mode prepared in the strong coherent state is not shown in the figure.) Under 
the condition that D has detected a state \G), the reduced state of the output 
signal mode becomes 

g[ = -TT2{Ugi®\F){F\U^fl), (29) 
P 

where 11 = |G) (G*!. The norm p = TriTr2 (...) is the probability of measuring the 
state \G) and thus the probability of generating the state g'l. Introducing the 
non-unitary (conditional) operator 

Y = {G\U\F) (30) 

acting in the signal-mode Hilbert space, we can rewrite ( pUD as 

g[ = -Yg,Y\ (31) 
P 

where the probability now reads p = Tti{Y giY"^). In practice, synchronized se- 
quences of light pulses could be fed into the input ports of the two-mode cou- 
pler, each pulse being prepared in the respective state, and the modes are thus 
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non-monochromatic ones. In what follows we assume that the U(l,l) and U(2) 
transformations do not vary with frequency within the spectral bandwidth of the 
pulses, so that the formulas given in Sec. |l| directly apply. 

In order to write F in a closed form, we first represent the states \F) and \G) 
in the form of 



\F) = F(St)|0) = ^i^Jt-|o), 

m=0 

oo 

\G) = G(6t)|0) = J]G„6t"|o) 



(32) 
(33) 



n=0 



(or with b"^ being replaced by c^) and substitute these expressions into (|30|). Using 
® and ([ 



applying the s-ordering rule ||T7 



mm m,n 



and introducing h = a^a, we obtain (see Appendix P) 
Ya = {0\G\b^)UAF{P)\0) 



(34) 



Yc 



\g^ {-R*T-^a) F{P*R*a)] T*-\PTY 

(0|Gt(ct)?7cF(ct)|0) 
= {g^ {-TZT*-^a^)F{Vnd))] {VTf, 

I ) sc 

where the respective ordering parameters are given by 



(35) 
(36) 

(37) 



Since |Tp<l<|Tp, we see that sa, sc > 1. Note that the arguments of F and 
G in (^) are the adjoints of the corresponding arguments in (^). 

Coherent displacements can be separated from the ordering procedure. To see 
this, we derive from ([19D and (pOf) together with ( p3l) and ( p4|) the transformation 
formulas for the displacement operators of the signal and idler mode. 





|T|2 + 1 




|r|2 + i 




|T|2 - 1 




|T|2 - 1 



Ua DM)b2{(3) Ui 

= bi[P(Ta - Rl3*)]b2[P\-Ra* + T/5)], 

lie bM)Hl3) ill 

= bi[V{Ta + 7^/3 )]b2[V{-TZ*a + T*/3)]. 



(38) 



(39) 
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Combining ( p8| ) and ( p9D with ( p5| ) and (p6D, respectively, yields 



f Pa* -Td*\ - f P*d* -T*a*\ 
\l = (li\G\c*)Dl(l3)Ucb.Aa)F(i*m 

i.e., a coherent displacement of the idler mode is equivalent to a corresponding 
coherent displacement of the signal mode. 

Since each trial in Fig. |I| yields a desired measurement outcome only with some 
probabilility, a lockable aperture S is needed in order to extract the properly 
transformed outgoing signal modes and dump the others. The desired output 
states are thus available at random times. It may however be demanded to 
provide them at certain times. In this quantum state memory M has 

to be used into which the pulses can be fed and released when desired. One 
possibility to realize M is offered by electromagnetically induced transparency 
. If, in particular, a pulse train of a certain repetition frequency is required. 



one may apply, e.g., an array of delays with variable optical path lengths which 
step by step adjust the waiting periods between the pulses to each other. Besides 
electromagnetically induced transparency, the cross-Kerr effect should also offer 
a way to realize a variable optical delay. In the latter case the outgoing signal 
mode in the state q'^ is mixed with a reference mode (index 3) prepared in a 
coherent state |a) at a cross-Kerr medium ^ = e''^"i"'3_ Yoi k — and |a| — ^cxd 
with = const., the reduced signal state becomes q'{ = TisCU g[\a) {a\tJ^) = 
giK|a| ni rai_ yjg^ choosiug |ap wc cau therefore control the refractive index 
and with it the optical length and time delay caused by the cross-Kerr medium. 



3 Preparation of single-mode quantum states 

3.1 Displaced photon adding for generating qubits 

In order to illustrate the general results derived in Sec. ^ let us first study the 
generation of a single qubit. When a parametric amplifier is fed with an idler 
pulse prepared in a single- mode coherent state \F) = \ai) and a single photon 
is detected, \G) = then (|35|) and ( ^OD yield the non-unitary (conditional) 
operator 

rj^-) = -RP*f*-^D{PR*-^al) 

x{PTy-^a^D{-T*R*-^al) (42) 
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(the index k = l is introduced for later purposes). Let us further assume that the 
signal input channel is unused, ^= |0)(0| (for notational convenience we omit the 
mode index). The outgoing signal pulse is then prepared in a state 



ff |0)(0|fft 



where 



1^) 



|0) + g|i) 



(43) 
(44) 



is a superposition of the vacuum and a single-photon Fock state. The parameter 
q = —PR/ai can be controlled by varying i? or ai. It is however convenient to 
choose \R\ and |ai| such that for a desired q the probability 



P 



(1) 



A 



(|i?|2+ |ain|T|-Vl"il 



(45) 



(111$) II = a/ ($1$)) of generating the qubit ( PD attains a maximum. This is the 
case for 



\R\ 





-2 + l)2 + 4|g 




-2 + 1) 


2 


Q 




2 





(46) 



The maximum values of p together with the corresponding values of |-R| and |tti| 
are shown in Fig. ^. 



1 

0.8 
. 6 
. 4 
0.2 



\R\ 







0.2 



. 4 



. 6 



0.8 



■ arctan |^| 



Figure 2: The maximized probability (^5|) of generating a qubit and the corresponding 
values of the parameters |i?| and |q;i| as functions of \q\. 



3.2 Repeated displaced photon adding for generating ar- 
bitrary superpositions of Fock states 

The scheme can be extended to the generation of an arbitrary superposition of 
Fock states 

N 

\^) = J2\n){n\^). (47) 

n=0 



Since the states of that type are completely determined by the N zeros of the 
Q-function, i.e., the N solutions Pi, - ■ ■ , Pn of the equation (\E'|/?) =0, they can be 
generated from the vacuum by alternate application of the coherent displacement 
operator and the creation operator, 



TV 



I*) 



k=l 

N 



n [oma^D^Pk)] |0), 

k=l 



(4J 



which may be realized repeating the procedure in Sec. ^]T] according to Fig. ^ 
The pulse prepared in the state is sent back through a ring resonator to 



Mo 




Ml 



I*) 



Figure 3: Scheme for preparing a travelling optical field in a quantum state (^). A 
parametric amplifier A is fed with a sequence of idler pulses prepared in appropriately chosen 
coherent states \ak)- They (and the pump pulses, not shown) arrive at A simultaneously with 
the produced signal pulse circulating in the ring resonator consisting of mirrors M2, M3 and 
a (removable) mirror Mi. The desired quantum state is generated if in each round trip of the 
signal field the detector D registers a single photon. 



the amplifier and used as signal input. Simultaneously, an idler pulse prepared 
in a coherent state \a2) is fed into the second input port of the amplifier. If 
the detector registers again a single photon, then the outgoing signal pulse is 
prepared in a state ~ Yj^^Yj^^\0). Under the condition that in each round trip 
a single photon is detected, after round trips the signal pulse is prepared in a 

|vi/)~rA|o) = ri^)---rf ri'^io). (49) 

If the preparation of the desired state has been successful, the mirror Mi in Fig. ^ 
can be removed in order to open the cavity and release the pulse. Otherwise the 
pulse is dumped in order to start the next trial from the very beginning, with the 
signal input port of the amplifier being unused. 



10 



Inserting (|i2| ) and rearranging the operator order such that the photon cre- 
ation operators are on the left of the exponential operators, we derive 



X 



TV 

n 

fc=i 



R 



l=k 



{pry 



|o) 



(50) 



(ttAT+i = 0; e'^ irrelevant phase factor). In the derivation of (0) we have used the 
relations D{a)f{a, a'^)D'^{a) = f{a — a, «"'' — «*) and a"'f{a, a'^)a~"' = f {a^^d, ad''). 
Comparing ( [50| ) with (^8]), we see that for 



N 



(3k 



(PT)^^Pat-T*a!^, 



or equivalently 



R* 



PR 



l=k 



{pry 



N 



Oik 



{PTY{PT) 



*N 



^|Tp'(A*-A!,i; 



(51) 



(52) 



l=k 



{Pn+1 = 0) the desired state is just realized. The probability of generating 
the state is given by 



P 



I^a|o)| 





R 


2N 


\{Nm\' 


2^ iV(Ar+3) 



exp 



N 

E 

k=l 



(53) 



as is seen by comparing the norms of both sides of ( pOD and inserting (^8]). It 
decreases rapidly with increasing N in general, so that the applicability of this 
method is effectively restricted to low numbers of round trips. 



3.3 Fock state generation 

In order to generate a Fock state \n), no feeding with idler modes is necessary. 



since all the (3k in (il) are zero. According to (|53D, the probability of detecting 
one photon at each of the n round trips is 



P 



^!|/2|2n|2^|-ri(7i+3) 



(54) 



and becomes maximal for fixed n if = 2/(n+ 1). Substituting this expression 
into (0) and applying Stirling's formula yields the asymptotic behaviour Pma.^ 
where a = \/27i/e and 6 = 2 — In 2, i.e., the probability of preparing 



ae 



-bn 



an n-photon Fock state decreases exponentially with increasing n. Note that 



11 



exactly the same asymptotic behaviour is observed when the state is generated 
by conditional measurement at a beam splitter array p|. 

However, the scheme considered here offers the specific feature of circumvent- 
ing the problem of low preparation probability. Since the idler field remains in 
the vacuum state, \ak) = |0), the circulating signal pulse is always in a Fock state 
whose number is simply the sum of detected idler photons. This means that every 
trial sooner or later results in the desired state \n), provided that the respective 
n is not skipped, and thus there is no need to wait for a sequence of n successive 
single-photon clicks. The idler detector is simply used for photon book-keeping 
and the cavity is opened in that moment when the sum of all detected idler 
photons has reached the desired value n. 



3.4 Influence of non-perfect photodetection and non-per- 
fect cavity feedback on Fock state generation 

If the idler detector and the feedback mirrors are not perfect, the situation be- 
comes more complicated. Let us assume that the state of the signal pulse after 
the jth round trip is a statistical mixture of Fock states described by a density 
operator g^^\ This pulse now enters the signal input port of the parametric am- 
plifier whose idler input port is unused, |F) = |0). When k outgoing idler photons 
are detected with efficiency ?7d, then the state g^^^' of the outgoing signal field is 
given by (E9]), where now 



Here, the symbol : : introduces normal ordering, and 



hkn{z) 



n 



\n—k 



(55) 



(56) 



|T|2p(A;) 



Inserting (|55|) into (|29|), we find that the prepared state is a mixture of Fock 
states with 

^ Y.h^^.i{r^,,)hU\T\-')Qf, (57) 
I 

where 

p{k) = |Tr2^6fc,„.K^D)&/„^(|Tr^)f?!/) (58) 

l,m 

is the probability of detecting k photons. Next, the pulse prepared in the state 
Q^^^' is fed back by the mirrors Mi,---,M3 into the signal input port of the 
amplifier. If the reflectances of the mirrors Mi,M2,M3 are given by 72-1, 7?.2, T^-s, 
the state of the pulse after the (j -|- l)th round trip is still a mixture of Fock states 
with 
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where r^p = |7^i7^2'7^3| determines the feedback efficiency. Combining ([59|) and 



we obtain the recursion relation 



\T\'^p{k 



'^■+1' l,n 

y<bk,^„i-n{VB)bni{\Tn0^£- (60) 

If we start from the vacuum state, Qnl = Sno, and measure the numbers kj of 
(outgoing) idler photons detected at the jth round trip, the circulating pulse 
therefore evolves into a mixture of Fock states, which depends according to (|60|) 
on the respective sequence {ki, k2, ■ ■ ■}. For this reason, some arbitrariness has 
to be introduced if the state evolution wants to be simulated. 

Let us therefore first consider the evolution of the mean photon number of 
the circulating pulse which can be obtained by considering the case when no 
measurement is performed in the idler output channel. ( |60D together with ?7d = 
and fcj+i = (and the initial condition gnl = 5no) then yields a thermal state 

'''' - wh-Mwhi) • 

where the mean photon number 

(n)(^)=,p|i^P ^^^'['^r-^ (62) 
r/FK I 1 

can be deduced from the recursion relation 

(n)(^'+i) = r]F [\TWny^^ + \R\^] . (63) 

If ?7f|^P > 1; then (n)*-^-* increases exponentially with A^, while for t^fIT^ < 1 a 
stationary value is observed, limAr^oo('^) =Vf\R\'^/ (1 — VfIT]"^). The critical value 
?7f|^P = 1 leads to a linear increase, (n)*^^^ = \R/T\'^N. 

To give an example, let us now consider the generation of the state |n = 4). 
Since the amplification is typically weak (|-Rp -C 1), a large number of round trips 
is likely and a high-quality cavity is demanded. Note that |-Rp = s(vac|d"'^a|vac)s 
= s(vac|6"l'6|vac)s characterizes the mean photon numbers of a two- mode squeezed 
vacuum |vac)s = f/A|0,0) that is generated by means of a parametric amplifier. 
In order to provide a rough estimate of the required efficiency rjp, we estimate 
the number N of cycles needed on average by identifying (n)^'^^ in ( |62D with n 
and assume a linear increase (n)^^^ = \R/T\'^N (because of the smallness of |-Rp), 
so that N \R\~'^n. We now estimate rj-p from the requirement that after that 
number of round trips in an empty resonator (i.e., without amplifier) an initially 
present photon can still be found with probability 1/2, thus rj-p = 2~l^l For n 
= 4 and = 3 ■ 10^^ this yields rjp ^ 0.999. We insert these quantities into (B^) 
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and (arbitrarily) consider the case when the first idler photon is detected if rip{n) 
exceeds 1, the second if rip{h) exceeds 2 and so on, until eventually the cavity 
is opened after detecting the 4th idler photon. The resulting density matrix is 
then obtained from (^). Examples are shown in Fig. ^ The plots confirm the 
sensitivity to cavity loss for small |-Rp. Note that the Mandel-parameter Q = 
{{Ahf)/{n) - 1 of the mixture plotted in (d) is (5 = -0.527. 
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Computer simulation of the preparation of a desired Fock state \n — A). 



The 



Figure 4: 

density matrix elements Qnn obtained according to (60) together with Qnn =Sno are shown for 
different feedback efficiencies yyp and detection efficiencies t^d, and = 3 • 10~'^. The cavity 
is assumed to be opened after detecting the 4th idler photon, which determines the number of 
round trips N . 

(a) r]F = 1 and ryo = 1 = 538). (b) r/p = 0.999 and 770 = 1 (A^ = 652). 
(c) r^F = 1 and 7/d = 0.7 (A^ = 636). (d) r^p = 0.999 and 770 = 0.7 (A^ = 788). 



4 Conclusion 

In this paper we have studied conditional quantum-state engineering at para- 
metric amplifiers and frequency converters, regarding each apparatus as being 
effectively a two-port device, whose action in conditional measurement can be 
described by a non-unitary operator Y defined in the Hilbert space of the sig- 
nal mode. We have presented Y for arbitrary quantum states of the incoming 
idler mode and arbitrary detected quantum states of the outgoing idler mode as 
s-ordered products of the operators that generate the quantum states from the 
vacuum, s being entirely determined by the devise parameters. 

To illustrate the results, we have proposed a scheme allowing the generation 
of arbitrary finite single-mode quantum states of travelling waves by a parametric 
amplifier equipped with a ring resonator as an optical feedback loop. We have 
applied the method to the problem of Fock state preparation, for which we have 
also addressed the influence of non-perfect photodetection and non-perfect cavity 
feedback. 
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A Derivation of Equation ( [2Sl ) 

The commutation relation between a and a''^ is invariant under a U(p,g) transfor- 
mation, i.e., for U^^ = G^U^G we have 

[(Ua)„(Ua)i] = [a,,at] = G,^. (64) 

After writing 

^giatHa ^ u-^Uae-'UtUt-iHU-Ua 

= u-i(Ua)e^(Ua)nu-tHu-){Ua)^ (55^ 

we choose U such that U~^''"HU~^ becomes diagonal, i.e., (U^^'''HU^^)a/x = 
(\J~^^'H.\J~^)x\Sx^. Applying the relation a'^f{a,a^)a~"' = f{a^^a,aa^), we see 
that the relation a^e'"'^>''*^ = e^'^'^A'*^e'*^^^"a^ holds. Inserting it into (^) with a 
= (U~^'''HU~^)aa, we obtain together with (|6^) 

^giatHa ^ u-i(Ua)ne'(""'"""^^^("'^^("'^^ 

A 

= U-ie^^'«V^(u-itHU-)(Ua) 

= giatHau-leiGU-^^HU-Ua 

= giatHagiGHp^^ (gg^ 

Note that e'*^^ is just a U(p,g) matrix, because e~'*^^ = G^e^'^^*^^G. 

B Derivation of Equation ( |S5l ) 

Applying (0), we first write 

6"f/A = t/A(?7l&^f/A)^" = UA{-P*Ra^ + P*ThY. (67) 
Using a"/(a, a''")a~"' = /(a^-'^d, ad"'"), we see from (^ that for |0) = |0)2 

(0|f/A = f*-\PTy-^'%Q\e^^\ (68) 
With the help of (|67D and (^) we now get 
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Applying 



J2 F^G*^T*'\PTy-^^^{0\e^^ 
^ FrnG*„T*-\PTy~^'^k\(j^^ Cf\ {P*Tf 

m,n,k \ / \ / 

X {R*T*-^ar-\-P*Rd^y~'' 

m,n,k ^ / \ / 

X {P*R*a)'^-^{-RT*'^a^Y~^{PTy-^ 
FmG*^{P*R*n~RT*-'y 



m,n,k 

X ■ 



X 



m\ in 
k)\k 



LL (Jj 



k 

with t = —1 (denoting antinormal order) and s = sa gives 



(69) 



J2 F^G*^{P*R*r{-RT*-^y 

m,n 

X {d^''d'^}s^f *-\PTy-^ 

J2 G*^{-R*T-^dy''Fm{P*R*d 

. m,n 

xf*-\PTy-^, 

from which (|35D directly follows. (|36|) can be derived analogously. 



(70) 
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